Context-free Grammars and Multivariate Stable 
Polynomials over Stirling Permutations 

William Y.C. Chen 1 , Robert X.J. Hao 2 and Harold R.L. Yang 3 

Center for Combinatorics, LPMC-TJKLC 
Nankai University, Tianjin 300071, P. R. China 

E-mail: ^hentQSnankai. edu.cn, 2 nalanxindao@163.com, 
3 yangruilong@mail. nankai. edu.cn 

o 

CN ■ Abstract 

Recently, Haglund and Visontai established the stability of the mul- 
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tivariate Eulerian polynomials as the generating polynomials of the Stir- 



"^^1 ling permutations, which serves as a unification of some results of Bona, 

Brenti, Janson, Kuba, and Panholzer concerning Stirling permutations. 
Let B n (x) be the generating polynomials of the descent statistic over 
Legendre-Stirling permutations, and let T n (x) = 2 n C n (x/2), where C„(x) 
' are the second-order Eulerian polynomials. Haglund and Visontai pro- 

, posed the problems of finding multivariate stable refinements of the poly- 

. nomials B n (x) and T n (x). We obtain context-free grammars leading 

to multivariate stable refinements of the polynomials B n (x) and T n (x). 
, Moreover, the grammars enable us to obtain combinatorial interpretations 

of the multivariate polynomials in terms of Legendre-Stirling permutations 
and marked Stirling permutations. Such stable multivariate polynomials 
provide solutions to two problems posed by Haglund and Visontai. 
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This paper presents an approach to the construction of stable combinatorial 
polynomials from the perspective of context-free grammars. The framework of 
using context-free grammars to generate combinatorial polynomials was pro- 
posed by Chen [5J. More specifically, we introduce the structure of marked 
1 Stirling permutations, and we find context-free grammars that lead to multivari- 

ate stable polynomials over marked Stirling permutations and Legendre-Stirling 
permutations. These multivariate stable polynomials provide solutions to two 
problems posed by Haglund and Visontai [13| in their study of multivariate 
stable refinements of the second-order Eulerian polynomials. 

Let us first review some backgrounds on the second-order Eulerian poly- 
nomials. These polynomials were first introduced by Gessel and Stanley [10] , 
which arc defined as the generating functions of the descent statistic over Stir- 
ling permutations. Recall that a Stirling permutation of order n is a permu- 
tation 7r = 7Ti7T2 • ■ ■ 7r2n-i7T2 n of the multiset {1, 1, 2, 2, . . . , n, n}, denoted by 
[n]2, which satisfies the following condition: if 7T; = ttj then irk > 7r^ whenever 
i < k < j. For 1 < i < 2n, we say that i is a descent of tt if i — In or iti > 7r$+i. 
Analogously, i is called an ascent of ir if i = 1 or 7Ti_i < 7T;. Let Q n denote the 
set of Stirling permutations of order n. Let C(n,k) be the number of Stirling 
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permutations of [n]2 with k descents, and let 

n 

C n (x) =Y J C{n,k)x k . 
fc=i 

Gessel and Stanley JTU] showed that 

±S(n + k,k)x" = ° n W , 

n=0 V ' 

where S(n,k), as usual, denotes the Stirling number of the second kind. The 
numbers C(n,k) are called the second-order Eulerian numbers by Graham, 
Knuth and Patashnik [11], and accordingly the polynomials C n (x) are called 
the second-order Eulerian polynomials by Haglund and Visontai |13j . 

The Stirling permutations were further studied by Bona [1] , Brenti [5] , Jan- 
son [14] and Janson, Kuba and Panholzer [15]. Bona p] introduced a statis- 
tic, called plateau, on Stirling permutations, and proved that ascents, descents 
and plateaux have the same distribution over Q n . Given a Stirling permuta- 
tion tt = 7T17T2 • • • 7T2„ 6 Q n , the index i is called a plateau of tt if tt{-\ = 7Tj. 
Analogous to that of the classical Eulerian polynomials, Bona Q] obtained the 
real-rootedness of the second-order Eulerian polynomials C n (x). 



Theorem 1.1 For any positive integer n, the roots of the polynomial C n (x) are 
all real, distinct, and non-positive. 

It should be noted that the real-rootedness of C n {x) is essentially the real 
rootcdness of the generating function of generalized Stirling permutations ob- 
tained by Brenti [5]. A permutation tt of the multiset {l ri , 2 1 " 2 , . . . , n r ™} is called 
a generalized Stirling permutation of rank n if tt satisfies the same condition as 
for a Stirling permutation. Let Q* denote the set of generalized Stirling per- 
mutations of rank n. In particular, if r% = = • ■ ■ = r n = r for some r, then 
tt is called an r-Stirling permutation of order n. Let Q n (r) denote the set of 
r-Stirling permutations of order n. It is clear that 1-Stirling permutations are 
ordinary permutations and 2- Stirling permutations are the Stirling permuta- 
tions. Brenti [5] showed that the descent generating polynomials over Q n have 
only real roots. 

Janson }14) defined the following trivariate generating function 
C n (x,y,z) = Y, x des My asc MzP lat M, 

where des(-7r), asc(7r), and plat(7r) denote the numbers of descents, the number 
of ascents, and the number of plateaux of n, respectively, and proved that 
C„(x,y, z) is symmetric in x, y, z. This implies the equidistribution of these 
three statistics derived by Bona. 

The symmetric property of C n (x,y, z) was further extended to r-Stirling 
permutations by Janson, Kuba and Panholzer [15] . For an r-Stirling permuta- 
tion, they introduced the notion of a j-plateau. For an r-Stirling permutation 
7r = 7Ti7T2 . . . 7T nr and an integer l^j'sCr — l,a number 1 ^ i < nr is called a 
j-plateau of w if 7T; = 7Tj+i and there arc j — 1 indices I < i such that 717 = m, i.e., 
the number 717 appears j times up to the i-th position of tt. Let j-plat(7r) denote 
the number of j-plateaux of tt. Meanwhile, define a descent and an ascent of tt 
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similar as ordinary permutations, and let des(7r) and asc(-7r) denote the number 
of descents and ascents of it. Janson, Kuba and Panholzer [Kj showed that the 
distribution of (des, 1-plat, 2-plat, . . ., (r — l)-plat, asc) is symmetric over the 
set of r-Stirling permutations. 

Based on the theory of multivariate stable polynomials recently developed 
by Borcea and Brandon [2HI], Haglund and Visontai [13] presented a unified 
approach to the stability of the generating functions of Stirling permutations 
and r-Stirling permutations. A polynomial /(z) G C[z] = C[z 1; z 2 , . . . , z m ] is 
said to be stable, if whenever the imaginary part lui(zi) > for all i then 
/(z) 7^ 0. Clearly, a univariate polynomial f(z) € K[z] has only real roots if 
and only if it is stable. 

For the case of univariate real polynomials, Polya and Schur |16j charac- 
terized all diagonal operators preserving stability or real-rootedness. Recently, 
Borcea and Branden [2H1] characterized all linear operators preserving stability 
of multivariate polynomials, see also the survey of Wagner [18]. This implies a 
characterization of linear operators preserving stability of univariate polynomi- 
als. 

A multivariate polynomial is called multiaffine if the degree of each vari- 
able is at most 1. Borcea and Branden showed that each of the operators 
preserving stability of multiaffine polynomials has a simple form. Using this 
property, Haglund and Visontai [T3] obtained a stable multiaffine refinement 
of the second-order Eulerian polynomial C n (x). Given a Stirling permutation 

7T = 7T17T2 ■ ■ ■ 7T 2 „ G Qn, kt 

A(w) = {i|7Ti_l < 7T,;}, 

D(w) = {i\TTi > 7T i+ l}, 

P(?r) = {i\-?Ti-l = 7Tj} 

denote the set of ascents, the set of descents and the set plateaux of it, respec- 
tively. Define 

C„(x,y,z)= 2J II x iri Yi II Z7r >- 

Haglund and Visontai [13] proved the stability of C„(x, y,z). 



Theorem 1.2 The polynomial C n (x, y, z) is stable. 



It is worth mentioning that, as observed by Haglund and Visontai, the recur- 
rence relation between C n _i(x, y, z) and C„(x, y, z) can be used to derive the 
symmetry of C„ (x, y, z), which implies the symmetry of C n (x, y, z) obtained by 
Janson, Kuba and Panholzer |15j . 

Moreover, Haglund and Visontai extended the stability of C„(x, y, z) to gen- 
erating polynomials of r-Stirling permutations by taking the j-plateau statistic 
into consideration. Let Pj(ir) denote the set of j-plateaux of n. Haglund and Vi- 
sontai |13| obtained the following multivariate stable polynomial over r-Stirling 
permutations 



£!„(x,y,zi,...,z r _i) 
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They also obtained a similar multivariate stable polynomial for generalized Stir- 
ling permutations. 

In view of the real-rootedness of C n {x) and its multivariate stable refinement 
C n (x, y,z), Haglund and Visontai posed the problem of finding multivariate 
stable polynomials as refinements of the generating polynomials of the descent 
statistic over Legendre-Stirling permutations. The Legendre- Stirling permuta- 
tions were introduced by Egge [9] as a generalization of Stirling permutations in 
the study of Legendre-Stirling numbers of the second kind. For any n > 1, let 
M n be the multiset {1, 1, 1, 2, 2, 2, • • ■ , n, n, n). A permutation n = Tri^ • • ■ ^3 n 
on M n is called a Legendre-Stirling permutation if whenever i < j < k and 
TTi = 7Tfc are both unbarred, then ttj > 7Tj. For a Legendre-Stirling permutation 
7r on M n , we say that i is a descent if either i = 3n or Hi > iti+i. Let -B n ,fc 
denote the number of Legendre-Stirling permutations of M n with k descents. 
Define 

2n-l 

B n {x) = B n,kX k . 

k=l 

Egge proved the real-rootedness of B n (x) . 

Theorem 1.3 For n > 1, B n (x) has distinct, real, non-positive roots. 

In order to derive a multivariate stable refinement of B n (x) , we introduce an 
approach of generating stable polynomials by a sequence of grammars. Based 
on the Stirling grammar given by Chen and Fu [7| , we find a sequence G\ , G%, . . . 
of context-free grammars to generate Legendre-Stirling permutations. We show 
that the formal derivative with respect to G n preserves stability by applying 
Borcea and Brandon's characterization of linear operators preserving stability. 
This leads to a multivariate stable refinement 5 n (x, y, z, u, v) of B n (x). On the 
other hand, according to the grammars, we obtain the following combinatorial 
interpretation 

£„(x,y,z,u,v) = ^2 n n y*i n ^ n ^ n u ^ 

7T iGX(n) iey(TT) ieZ(n) i£U(Tr) i£_V(ir) 

The real-rootedness of B n (x) is a consequence of the stability of £> n (x, y, z, u, v) 
by setting Vi = yi = y and Xi = Zi — Ui — 1 for < i < n. 

Haglund and Visontai also posed the problem of finding multivariate stable 
refinements of the polynomials T n (x), which are given by 

T n (x) = 2 n C n (|) = £ 2 n ' k C(n, k)x k , (1.1) 

k 

where C(n, k) and C n (x), as before, denote the second-order Eulerian numbers 
and the second-order Eulerian polynomials respectively. The polynomials T n (x) 
were introduced by Riordan [T7] . 

In light of the relation (jl.ip between T n (x) and C n (x), we introduce the 
structure of marked Stirling permutations and the following multivariate poly- 
nomials 

7 i n(x,y,z)=^ Yi n y7Ti n 

t ieD(ir) i£A(iv) ieP(Tr) 

where n ranges over marked Stirling permutations of [n}2- We shall show that 
the polynomials T„(x, y,z) are stable. The polynomial T n (x) becomes the spe- 
cialization of T„(x, y, z) by setting Xi = Zi = 1 and yi = x for < i < n. This 
implies that T n (x) is real- rooted. 
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This paper is organized as follows. In Section [21 we give an overview of dif- 
ferential operators associated with context-free grammars. We find context-free 
grammars to generate the polynomials C n (x, y, z). In Section [31 we obtain 
context-free grammars that lead to the multivariate generating polynomials 
-B„(x, y, z, u, v). In Sectional we introduce the structure of marked Stirling 
permutations, and we give context-free grammars to generate the multivariate 
polynomials T„(x,y, z). In Section [5l based on Borcea and Branden's charac- 
terization of stability preserving linear operators, we present an approach to 
proving the stability of polynomials generated by context-free grammars. In 
particular, we prove the stability of multivariate polynomials i?„(x, y, z, u, v) 
and T„(x,y, z). 

2 Context-free grammars 

In this section, we give an overview of the idea of using context-free grammars 
G to generate combinatorial polynomials and combinatorial structures as de- 
veloped by Chen [BJ. A context-free grammar G over an alphabet A is defined 
to be a set of production rules. Roughly speaking, a production rule means to 
substitute a letter in the alphabet A by a polynomial in A over a field. Given a 
context-free grammar, one may define a formal derivative D as a linear operator 
on polynomials in A, where the action of D on a letter is defined by the substi- 
tution rule of the grammar and the action of D on a product of two polynomials 
u and v is defined by the Leibnitz rule, that is, 

D(uv) = D(u)v + uD(v). 

Many combinatorial polynomials can be generated by context-free gram- 
mars. Meanwhile, context-free grammars can be used to generate combinatorial 
structures. More precisely, one may use a word on an alphabet to label a com- 
binatorial structure such that the context-free grammar serves as the procedure 
to recursively generate the combinatorial structures. Such a labeling of a com- 
binatorial structure is called a grammatical labeling in [7]. 

For example, the grammar 

G = {a -> ab, b — > b} 
is used in [6J to generate the set of partitions of [n] and the Stirling polynomials, 

n 
i=0 

where S(n, k) denotes the Stirling number of the second kind. For a partition P, 
we label a block of P by letter b and label the partition itself by letter a, and we 
define the weight of a partition by the product of its labels. So a partition P with 
k blocks has the weight w(P) = ab k . For example, the partition {{1, 2}, {3}} is 
labeled as follows 

{1,2}{3} 

b b a . 

In the above notation, we write a partition P = {Pi, P2, . . . , Pk} of [n] in such 
a way that the blocks are ordered in the increasing order of their minimum 
elements. Moreover, we put the letter a at the end of the partition. 



Using the above grammatical labeling of a partition, we deduce that 

n 

D n (a) =J2 w ( p ) =J2 S ( n > k ) abk ' t 2 ' 1 ) 

P k=l 

Many properties of the Stirling polynomials follow from the above expression in 
terms of the differential operator D with respect to the grammar G. 

Let us explain how the grammar works for the generation of partitions. For 
n = 1, there is one partition of [1], that is, {{1}}, whose label is ab. Assume 
that we have generated all the partitions of [n — 1] by applying the operator 
D n ~ 2 to {{1}} with the initial grammatical labeling. 

Let us give an example to demonstrate the action of the differential operator 
D with respect to the grammar G to a partition of [n] with the aforementioned 
grammatical labeling. Consider the following partition of {1, 2, 3, 4, 5, 6} 

{1,3,6}{2,5}{4} 

b b b a . 



If we apply the substitution rule to the letter a, then we get ab which we 
rewrite as 6a, where a still serves as the label of the new partition, and b stands 
for a new block {7}. In this case, we get a partition 

{1,3,6}{2,5}{4}{7} 

b b b b a . 



If we apply the substitution rule to the second letter b, then we get b. In 
this case, we insert the element 7 in the second block, and we are led to the 
following partition with consistent grammatical labeling 

{1,3,6}{2,5,7}{4} 

b b b a . 



Starting with the empty partition with label a, we get 

{i} 

D(a) = b a, 

{1}{2} {1,2} 

-D (a) = b b a + b a, 

{1}{2}{3} {1}{2,3} {1,3}{2} {1,2}{3} {1,2,3} 

Zr(a) = b b ba+ b b a+ b ba+ b ba+ b a. 



Without considering the combinatorial structures during the applications 
of the differential operator D, we may directly compute D n (x) to derive the 
Stirling polynomials S n {x). 

As the second example, we consider the context-free grammar 

G = {x -> xy, y -> xy} 

introduced by Dumont [8] which is used to compute the Eulerian polynomials. 
For a permutation tt = 1:1^2 ■ ■ ■ i"n of [n], let 

A(n) = {i I 7Ti_l < TTi}, 

D(ir) = {i\n> 
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denote the set of ascents and the set of descents of it, respectively. Here we 
set ttq = 7r n+ i = 0. In other words, for any permutation it of [n], 1 is always 
an ascent and n is always a descent. An element Tii is called a descent top of 
7r if i <E D(ir), and iri is called an ascent top if i € A(tt), see Haglund and 
Visontai (T3] . 

The grammatical labeling of a permutation tt is defined as follows. If 7r, is 
an ascent top of 7T, then we label with the letter x. If 7ri is a descent top, 
then we label 7Tj by the letter y. For this labeling, the weight of 7r is given by 

Then for n ^ 1, we have 

D n (x) = £ «;(*) 

where A(n, m) is the Eulerian number, namely, the number of permutations of 
[n] with m descents, see Dumont [3]. 

For n = 1, there is only one permutation of [1], that is 1, whose label is xy. 
Assume that we have generated all the permutations of [n — 1] by applying the 
operator D n ~ 2 to 1. 

Next we give an example to illustrate the action of I? on a permutation of 
[6] . Take a permutation 

3 2 5 6 4 1 

x y x x y y y . 

If we apply the substitution rule x — > xy to the third letter x, we insert 7 after 
5. As for the grammatical labeling, we keep all the labels and assign the element 
7 a new label y as if it comes from the substitution rule x — > xy. Indeed, it is 
easily checked that what we get is a permutation with a consistent grammatical 
labeling, namely, 

3 2 5 7 6 4 1 

xyxxyyyy. 

Similarly, if we apply the substitution rule y —> xy to the second letter y, 
then we insert 7 after 6. In this case, we need to change the label of 6 from y 
into x, and assign y to the new element 7. In other words, the label y becomes 
xy just like the substitution rule. So we get the following permutation with a 
grammatical labeling, 

3 2 5 6 7 4 1 

xyxxxyyy. 

Indeed, the above examples indicate that permutations of [n] and the Eu- 
lerian polynomials A n (x) can be generated by the operator D associated with 
the grammar G. 

In order to generate combinatorial structures with more parameters, we may 
use a sequence of grammars. Let us consider the the multivariate refinement of 
Eulerian polynomials ^4„(x,y) introduced by Haglund and Visontai |13j . which 
involve the sets of ascent tops and descent tops, not just the numbers of ascents 
and descents. More precisely, 

A»(x,y)= II x *' II y7r '- 

7rSS„ i£A{-n) iGD(-!r) 

We shall introduce a sequence of grammars {G n } to generate the multivariate 
polynomials A„(x,y). 



n 

= A(n,m)y m x n+1 - 

m—1 



7 



For n > 1, define 

y n , < i < n}, 

and denote by D n the formal differential operator with respect to G n . The mul- 
tivariate polynomials A n (x, y) can be generated by the sequence of grammars 
G n . 

Theorem 2.1 For n > 1, we have 

D n D„- 1 ■ ■■D 1 (x ) = A„(x,y). 

Proof. We define the grammatical labeling of a permutation n as follows. For 
a permutation it, if iri is an ascent top, we label 7Ti_i by the letter x Vi ; if TTi is 
a descent top, we label 7Tj by the letter y 7ri . So the weight of 7r is given by 

w(tt) = ii, J] y-Ki- 

iGA(-ir) i£D(7r) 

We proceed to show by induction that D n D n -\ ■ ■ ■ Di(xq) equals the sum 
of the weights of permutations of [n]. For n — 1, the theorem is valid since the 
weight of the permutation 1 is x\y±. Assume that the theorem holds for n— 1, 
that is, 

D n -x ■ ■ ■ Di(x ) = ^2 w(tt). 

We now use an example to illustrate the action of D n on a permutation of [n— 1]. 
Let 7T = 325641. The grammatical labeling is as follows 

3 2 5 6 4 1 

^3 2/3 x 5 x 6 y 6 y A y x . 

If we apply the substitution rule xq — > x 7 y 7 to the letter xq, we define the 
action as the insertion of 7 immediately after 5. The labels of 5 and 7 will be 
changed to x 7 and 2/7 as given by the grammar. It is not hard to see that the 
permutation we obtain has a consistent grammatical labeling, 

3 2 5 7 6 4 1 

£3 2/3 ^5 x 7 y 7 y 6 y 4 yi . 

Similarly, if we apply the substitution rule y$ — > 272/7 to the letter 2/6, we obtain 
a permutation with a consistent grammatical labeling 

3 2 5 6 7 4 1 
£3 2/3 ^5 X 6 X 7 J/7 2/4 2/1 • 

It is clear that all permutations of [n] can be obtained this way. So we conclude 
that 

£>„£>„-! • • • Si (so) =Dnl J2 W M = 

Hence the theorem holds for all positive numbers n by induction. I 

For 71 = 0, the empty permutation is labeled by xq. The values of A n (x, y) 
for n = 1,2,3 are given below. 

Dxixo) = x x 2/1, 

„ _. , . 21 12 

D 2 D 1 (x ) = x 2 2/2 2/1 + x\ x 2 2/2, 

321 231 213 312 

D 3 D 2 Di(x ) = x 3 2/3 2/2 2/1 + x 2 x 3 2/3 2/1 + x 2 2/2 x 3 2/3 + 2:3 2/3 x 2 y 2 

13 2 12 3 

+ xi x 3 2/3 2/2 + xi x 2 x 3 2/3 ■ 
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Let us now consider the grammar to generate Stirling permutations. Chen 
and Fu [7] showed that the grammar 

G = {x -> x 2 y,y^r x 2 y} 

can be used to generate Stirling permutations. Let D denote the differential 
operator associated with the grammar G. It has been shown in [7] that 

n 

D n (x) =x^2 C(n, m)x 2n - m y m , 

where C(n, m) denotes the second-order Eulerian number. Notice that 

D n (x)\ x=1 = C n (y), 

where C n {y) is the second-order Eulerian polynomial. 

The grammatical labeling of a Stirling permutation is defined as follows. For 
a Stirling permutation tt, if i € D(tt), we label tti by y; if i £ A(tt) or i £ P(tt), 
we label 7r 2 ;_i by x. For example, the Stirling permutation 7r = 233211 has the 
following grammatical labeling 

2 3 3 2 1 1 

x x x y y x y . 

Next we show that one can use a refinement of the grammar G to derive the 
multivariate polynomials C n (x, y, z) of Haglund and Visontai |13j . Recall that 

c„(x, y ,z)= n n y ^ n 

tt6Q„ ieA(?r) iED(n) ieP(7r) 

As a refinement of the grammar G, we define 

G n {.t ? Y x n y n z n: y.[ y x n y n z n , Zj ^ x n y n z n: ^ z <C 7^}. 
and we denote by D„ the differential operator associated with the grammar G n . 

Theorem 2.2 For n ^ 1, we have 

D n D n -i ■ ■■Dx{zq) = C„(x,y,z). 

Proof. First, let us define the grammatical labeling of a Stirling permutation n. 
For a Stirling permutation tt, if ^ is an ascent top, we label 7Tj_i by the letter 
x 1Ti ; if TTi is a descent top, we label iri by the letter y 7ri ; and if iTi is a plateau, 
we label iii-i by the letter z„ t . For this labeling, the weight of tt is given by 

w(ir) = x 7Ti y^ z Vi . 

i£A(ir) i£D(-k) ieP(-ir) 

We aim to show that D n D n -i ■ ■ ■ Di(zo) equals the sum of weights of Stirling 
permutations of [71)2- Let us use induction on n. The theorem is obvious for 
n = since the weight of the empty permutation is Zq. Assume that the theorem 
holds for n — 1, that is, 

D„_i • • •D 1 {z a ) = ^ w ( 7r )- 

7T£Q„_1 
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Let us use an example to demonstrate the action of D on a Stirling per- 
mutation of [n — 1] 2 . Let 7r = 233211. The grammatical labeling of tt is as 
follows 

233211 
x 2 x 3 z 3 y 3 y 2 zi yi . 

In general, if we apply a substitution rule of G4 to any letter in 7r, wc get £41/4214. 
Here we insert the two elements 44 after the element whose label is replaced by 
the substitution rule, and we use the labels X4, 2/4 and 2:4 to relabel the three 
elements that are affected by the substitution. For example, if we apply the 
substitution rule x 2 — > x^y^Zi to the above Stirling permutation, then we get a 
Stirling permutation with the following grammatical labeling 

44233211 

24 z A y 4 x 3 z 3 y 3 y 2 zi y 1 . 

It is easily seen that the application of any substitution rule of G n to any Stirling 
permutation of [n— l] 2 leads to a Stirling permutation of [n] 2 with a consistent 
grammatical labeling. Hence we deduce that 

D n D n -i---D 1 (z )=D n l »W =E W W- 

Thus, the theorem holds for n. This completes the proof. I 
For n = 0, the empty permutation is labeled by z . The values of the 
polynomials C„(x, y,z) for n = 1,2 are as follows, 

Di(zo) = xx z\ yi, 

2211 1221 1122 

D2D 1 {z ) = x 2 z 2 y 2 zi y\ + x\ x 2 z 2 y 2 y-y + xi zi x 2 z 2 yi ■ 

We shall give further refinements of the above two sequences of grammars 
as solutions to the problems of Haglund and Visontai [13] . On one hand, wc 
use these refined grammars to construct multivariate polynomials for Legcndre- 
Stirling permutations and marked Stirling permutations. On the other hand, 
we use the grammars to construct stability preserving operators leading to the 
stability of the multivariate polynomials. 

3 Legendre-Stirling permutations 

In this section, we introduce several statistics on Legendre-Stirling permutations 
of 

M n = {1,1, 1,2,2,2,- •• ,n,n,n}. 

In terms of these statistics, we obtain multivariate polynomials B n (pc, y, z, u, v) 
as refinements of B n (x). In fact, the combinatorial construction of the multi- 
variate polynomials -B„(x, y, z, u, v) is obtained from further refinements of the 
grammars to generate permutations and Stirling permutations with respect to 
the numbers of descents. Using these grammars, we derive the combinatorial 
interpretation by giving a suitable grammatical labeling. In Section [5l we shall 
use grammars to prove the stability of B n (x, y, z, u, v). This leads to a solution 
to the problem of Haglund and Visontai. 
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Let L n denote the set of Legendre-Stirling permutations of M n . For a 
Legendre-Stirling permutation tt G L n , define 



X(n) 


= {* 


7Tj 


-1 < 7Tj 


7Tj is unbarred and appears the first time}, 


Y(n) 


= 


7T, 


> + 1 


and Tti is unbarred}, 


Z{tt) 


= {* 




-1 ^ 7Tj 


7Tj is unbarred and appears the second time} 


U(tt) 


= 


7T, 


-1 < 7T, 


and iTi is barred}, 


V(tt) 


= {* 




> 7Ti+l 


and TTi is barred}. 



Here we set ttq = 7T3„+i = 0. 

For example, let tt = 112233231. Then we have X(ir) = {2,4,5}, Y(n) = 
{6,9}, Z(tt) = {6}, [/(tt) = {1,3,8} and V(tt) = {8}. 

Define 

5„(x,y, z,u,v) = ^ u Ti v m . 

7r£L„ ieX(Tr) ieY(ir) iEZ(n) ieU{n) i£V (tt) 

We define the grammars {G n } as follows, 

G 2n -i = {xi,yi,Zi,Ui,Vi -> u n v n ,0 s^i < n}, 
G 2n = {xi,Vi, Zi, Ui, Vi -t x n y n z n , < i < n; 

u n ^ x n z n u n ,v n y x n y n z n }. 

Notice that G2n-i is a refinement of the grammar 

G = {x -> xy,y -> xy}. 
and is a refinement of the grammar 

G = {x ->■ x 2 y,y ^ x 2 y} 

The grammatical labeling of a Legendre-Stirling permutation is defined as 
follows. Let tt be a Legendre-Stirling permutation on M n . For i G ^(7r), 
i G Z(tt) or i G U(tt), we label 7r 2 ;_i by the letter respectively; 
for i G ^(tt) or i G V(7r), we label i\i by the letter j/,^ or v Wi , respectively. 
For example, the above Legendre-Stirling permutation tt = 112233231 has the 
following grammatical labeling 

112233231 

Mi xi u 2 x 2 x 3 z 3 y 3 u 3 v 3 yi . 

The following theorem shows that the polynomials B n (x, y, z, u, v) can be 
generated by the grammars G„ . 

Theorem 3.1 For n ^ 1, we have 

D 2 nD 2 n-i ■ ■ ■ D^xq) = B„(x, y, z, u, v). (3.1) 
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Proof. We use induction on n. The case for n = is obvious since the empty 
permutation is labeled by xq. Assume that the theorem holds for n — 1, that is, 

D 2n -2 ■■D 1 (x ) = W W- ( 3 - 2 ) 

7r6L„-i 

Note that any Legcndre-Stirling permutation of M n can be obtained from 
a Legendrc-Stirling permutation of M n _i through two operations: (1) Insert a 
barred element n; (2) Insert two elements nn. We use an example to show that 
the operators Z?2n-i and Dm correspond to these two operations. 

Consider the Legendrc-Stirling permutation 7r = 112233231, whose gram- 
matical labeling is given by 

112233231 

ui xi u 2 x 2 x 3 z 3 y 3 u 3 v 3 yi . 

The first operation is just the procedure of generating permutations. In 
general, if we apply a substitution rule of GV to it, we always get U4V4. Here we 
insert 4 after the element whose label is replaced by the substitution rule. At 
the same time, we relabel the two involved elements by the letters 1*4 and W4. 
For example, if we apply the substitution rule z 3 — > U4V4 to tt, then we obtain 
a Legendre- Stirling permutation with a consistent grammatical labeling 

I122343231 
ui x\ u 2 x 2 x 3 Ui Vi y 3 u 3 v 3 y\ . 

As for the second operation, consider the above Legendre-Stirling permuta- 
tion a = 1122343231. The two substitution rules M4 — > x^ZiUi and V4 — > x^y^Zi 
of Gg correspond to the operations of inserting two elements 44 before 4 or after 
4, respectively. So we get two Legendre-Stirling permutations 112234443231 or 
112234443231. 

Next we consider the rest of substitution rules of Gg. If we apply any of 
the remaining substitution rules of G$ to a, we always get CC42/4Z4. Here we 
insert two elements 44 into a between 7^ and 7r J+ i, which are elements less than 
4. For example, by applying the production rule 112 — > x^y^Zi, we obtain the 
Legendre-Stirling permutation 

1 14422343231 

Wi Xi x 4 z 4 2/4 x 2 x 3 u A w 4 y 3 u 3 v 3 yi . 

It can be checked that any of applications of the substitution rules of Gs to a 
leads to consistent grammatical labelings. Moreover, it can be verified that the 
action of D2 n D2 n -i 011 the Legendre-Stirling permutations of M„_i generates 
all the Legendre-Stirling permutations of M n . So we conclude that 

D2 n D 2n -i ■ ■ ■ Dt(xo) = 2J w(tt). 

Then the theorem follows by induction. I 

For n = 0, the empty permutation is labeled by xq, and i?i(x, y, z, u, v) is 
calculated as follows, 

Di(xq) = Ml Vx, 
_. , , 111 111 

D 2 Dx{xo) = x\ zx ux vx + ux xx z Y yx ■ 
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4 Marked Stirling permutations 



In this section, we introduce the structure of marked Stirling permutations, 
and we define several statistics in order to construct multivariate polynomi- 
als T„(x, y,z) as refinements of T n {x). We also give a sequence of grammars 
to generate T„(x, y,z) as well as marked Stirling permutations with suitable 
grammatical labelings. By using the grammars, the stability of T n (x, y,z) can 
be established in Section O This gives a solution to the problem of Haglund 
and Visontai concerning a stable refinement of T n (x) . 

A marked Stirling permutation is defined by the following marking rule. 
Given a Stirling permutation tt = tt\tti ■ ■ ■ 7T2„, if TTi is an element of tt such that 
ni occurs the second time in tt and 7Tj < 7r i+1 , then we may mark the element 
7Tj. We denote a marked clement i by i. A marked Stirling permutation is a 
Stirling permutation with some elements marked according to the above rule. 

For example, there is only one marked Stirling permutation of [1)2: 11, 
whereas there are four marked Stirling permutations of [2] 2: 

2211,1221,1122,1122. 

Let Q n denote the set of marked Stirling permutations of \n\2- We use 
A(tt), D(tt), P(tt) to denote the set of descents, the set of ascents and the set 
of plateaux of tt. More precisely, given a marked Stirling permutation tt = 

7Ti7r 2 • • • 7T 2n £ Q n , let 

A{tt) = {i I 7Tj_i < TTi}, 

D(tt) = {i I TTi > 7T i+ l}, 

P(ir) = {i I 7Ti_i = 7Tj} 

denote the set of ascents, the set of descents and the set of plateaux of tt, 
respectively. Let T(n, m) be the number of marked Stirling permutations of 
[n\2 with m descents. It follows from relation (jl.ll) that 

n 

T n (x) = T(n,m)x m . 

Note that Riordan [17] introduced the polynomials T n (x) and proved that 
T n {l) equals the Schroder number, namely, the number of series-reduced rooted 
trees with n + 1 labeled leaves. 

We shall prove that the polynomials T n (x) can be generated by the grammar 
G defined by 

G = {x^x 2 2 /,y^2x 2 2 /}. 

The grammatical labeling of a marked Stirling permutation can be described as 
follows. Let tt be a marked Stirling permutation of [n]2- If i 6 D(tt), we label 
TTi by y. If i £ A(tt) or i 6 P{k), we label tti-i by x. The weight of a marked 
Stirling permutation tt of [n]2 with m descents is given by 

w(tt) = x 2n+1 - m y m . 
Theorem 4.1 For n 1, we have 

n 

D n (x) = T(n,m)x 2n - m+1 y m . 

m— 1 
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Setting x = 1, we have 

D n (x)\ x=1 =T n (y). 

Proof. We aim to show that D n (x) equals the sum of the weights of marked 
Stirling permutations of [n]2- We use induction on n. The case for n = 
follows from the fact that the weight of the empty permutation is x. Assume 
that the theorem holds for n — 1, that is, 

D n ~\x)= HO- 

We now use an example to demonstrate the action of D on a marked Stirling 
permutation of [n— 1]2- Let tt = 122331 with the following grammatical labeling 

1 2 2 3 3 1 

x x x x x y y . 

If we apply the substitution rule x — > x 2 y to the fourth letter x, then we 
insert the two elements 44 after 2. We keep all the labels except that we assign 
the labels x and y to the two new letters 44. It is not difficult to see that the 
generated marked Stirling permutation has a consistent grammatical labeling 

12244331 

xxxxxyxyy. 

If we apply the substitution rule y — > 2x 2 y to the first letter y, then we 
insert 44 after the second element 3. We change the label of the second element 
3 from y into x and assign x and y to the two new elements 44. According to 
the marking rule, the second element 3 may be marked or unmarked. These 
two choices correspond the coefficient 2 in the substitution rule y — > 2x 2 y. So 
we are led to the following two marked Stirling permutations with consistent 
grammatical labelings, 

12233441 

xxxxxxxyy, 

and 

12233441 

xxxxxxxyy. 

It can be verified that the above process generates all marked Stirling permu- 
tations of [n]2- It follows that 

D n (x)=D(D n ~ 1 (x))=D I ]T w(n) j = £ w(a). 

\7reQ„_i / cr£Q„ 

Hence the proof is complete by induction. I 

As a multivariate refinement of T n (x), we define the following generating 
polynomial of marked Stirling permutations of [n]2, 

r„(x,y,z)= ^2 n n y** n 

Let 

G n {.t^ Y x n y n z n: y^ y 1x n y n z n: y x n y n z n7 ^ i <C tt-}. 

The grammatical labeling of a marked Stirling permutation can be described as 
follows. For a marked Stirling permutation n of [n]2, if i € we label 7Tj_i 
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by Xi\ if i G D(ir), we label 71".; by y^ and if i £ f (vr), we label 7Ti_i by Zj. Then 
the weight of n equals 

w(ir) = x 7Ti y Xi z Vi . 

i£A(tt) ieD(ir) ieP(ir) 

The following theorem shows that the polynomials T„(x, y, z) can be gener- 
ated by the grammars G n . 

Theorem 4.2 For n > 1, we have 

£>„£>„_!•■•£>! (zo) = T„(x,y,z). 

The proof of the above theorem is analogous to that of Theorem 14. II Hence 
the details arc omitted. Here we use an example to illustrate the action of D 4 
to the above marked Stirling permutation n = 122331 with the grammatical 
labeling 

1 2 2 3 3 1 

X\ x 2 z 2 x 3 z 3 y 3 yi . 

If we apply the substitution rule 23 — > X4J/4Z4 of G4 to the letter x 3 , then we 
insert the two elements 44 after 2 to get a marked Stirling permutation with 
the following consistent grammatical labeling 

12244331 

xi x 2 z 2 x± z± y A z 3 y 3 yi . 

If we apply the substitution rule y 3 — > 2X4J/4Z4 of G4 to the letter y 3 , then we 
insert 44 after the second element 3 to get the following two marked Stirling 
permutations with consistent grammatical labelings 

12233441 

Xi x 2 z 2 x 3 z 3 X4 z 4 y 4 j/i, 



and 

12233441 

x x x 2 z 2 x 3 z 3 x 4 z 4 y 4 J/! . 

For n = 0, the empty permutation is labeled by zq. For n = 1, 2, T^(x, y, z) 
are given below, 

£>i(z ) = xx Zx yx, 

. . 2211 1221 1122 

D 2 Dx{z ) = x 2 z 2 2/2 zx yx + xi x 2 z 2 y 2 yx + x x Zx x 2 z 2 y 2 

1122 
+ Xx Zx x 2 z 2 2/2 • 

5 Grammars preserving stability 

In this section, we prove the stability of the multivariate polynomials i?„(x, y, z, u, 
and T„(x, y, z) based on context-free grammars and the characterization of sta- 
bility preserving linear operators due to Borcea and Branden [3]. 

Our idea of proving the stability of the polynomials by a sequence of context- 
free grammars {G n } goes as follows. Since the initial polynomial x is stable, if 
D 1} D 2 , . . . , D n , . . . preserve stability, then D n D n ^x ■ ■ ■ Dx(x) is stable. If D n is 
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not stability preserving, then we try to find a sequence of stability preserving 
operator {T n } such that 

T„T„_i . ..Ti(x) = D n D n _i . ..Di(x). 

If such operators T n exist, then we reach the conclusion that the multivariate 
polynomials D n D n -i ■ ■ ■ D\(x) are stable. 

Note that i?„(x, y, z, u, v) and T„(x, y, z) are all multiaffine polynomials in 
the sense that the degree in each variable is at most 1. In order to construct 
the stability preserving operators T n based on the grammars G„, we consider 
some equivalent forms of production rules when we restrict our attention to 
multiaffme polynomials. 

For example, let 

G n = {a — > ab n , b{ — > b ni < i < n}. 

Observe that as far as the computation is concerned, the formal differential 
operator D n with respect to G n is in accordance with the following operator 

n 

T n = b n (l + J2d/d h ), 
i=i 

when they are applied to certain polynomials. Thus we obtain 

T„T„_i . . .Ti(a) = D n D n ^ 1 . . . Di(a). 
However, D n and T„ are different operator in general, since 

D n (a + b 1 )^T n (a + b 1 ). 

For multiaffme polynomials, the characterization of stability preserving op- 
erators is simpler than that for the general case. For the purpose of this pa- 
per, we only need the following sufficient condition to prove the stability of 
-B n (x, y, z, u, v) and T n (x, y,z), see Borcca and Branden [3]. 

Lemma 5.1 Let f £ C[zi, Z2, ■ ■ . , z„] be a stable multiaffine polynomial and let 
T denote a linear operator acting on the polynomials in C[zi, Z2 1 ■ ■ . , z n ]. If 

T {X\( Zl + w ')^ e C[zx, z 2 , ■ ■ ■ , z n ,wi, . . . ,w n ] 

is stable, then T(f) is either stable or identically 0. 

Next we show how to prove the stability of polynomials generated by context- 
free grammars. Let us consider the multiaffine polynomials C„(x, y,z) defined 
by Haglund and Visontai [13]. Let 

and let G n denote the differential operator associated with the grammar G n . 
Let f n = D n D n -i ■ ■ ■ Di(zq). From the grammatical labelings, it is clear that 
/„ is multiaffine. We wish to prove the stability of /„ by induction on n. Since 
z is stable, it suffices to prove that the operator D n+ i preserves stability of 
multiaffine polynomials. 
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Let 

n 

F = Y[{xi + Wi)(yi + Vi)(zi + Ui). 
By Lemma it suffices to check the stability of D n+ i(F), that is, 



D n+ i(F) = x n+ iy n+1 z n+1 F^ 

i=0 



111 



y, 



is stable. 

If Xi, yi, Zi, Wi, Vi and Ui have positive imaginary parts for all < i < n, 
then 



=o 



Wi y t 



has negative imaginary part. Thus, 

D n+ i(F) = x n+ iy n+ iz n +iF£ / 0. 

Hence D n+1 (F) is stable. By Lemma [5.11 we find that D n+1 (f n ) is a stable 
polynomial. So we conclude that 

fn+i = F) n+ iD n D n ^i ■ ■ ■ Di(zo) 

is stable. 

The stability of A n (x, y) can be proved in the same way. Indeed, let 

y n , ^ i < n}, 

and let D n denote the differential operator with respect to G n . It turns out 
that the operator D n preserves the stability of multiaffinc polynomials. 

It is worth mentioning that the formal differential operators used in the above 
two examples arc essentially equivalent to the operators given by Haglund and 
Visontai [13] in their proofs of the stability of C„(x, y,z) and A n (x,y). 

Next we construct stable multivariate refinements of S n (x), the Stirling poly- 
nomials. Recall that the grammar 

G = {a ->• ab, b — > b} 

generates the polynomials S n (x). Define 

G n = {a ab n ,bi -> b n , 1 < i < n}, 

and let D n denote the formal differential operator associated with G n . We 
define the grammatical labeling of a partition as follows. For a partition P = 
{Pi,P 2 , . . . , Pfc}, we label the partition itself by the letter a and label a block 
Pi by the letter b m , where m is the maximum element in Pj. Then the weight 
of P is given by 

fe 

w{P)=aY[b mt , 

i=i 

where m, is the maximum element in Pj. Denote by S n (a, b) the sum of weights 
of partitions of [n]. The next theorem shows that S n (a, b) can be generated by 
G n - However, in this case, the differential operator D n associated with G n is 
not stability preserving even for multiaffinc polynomials. Instead, we shall find 
an equivalent operator T n that preserves stability for multiaffine polynomials. 
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Theorem 5.2 For n > 1, we have 

D n D n _ 1 .,.D 1 (a) = S n (a,b). 

The proof of the above theorem is analogous to that of (|2.ip . Here we use the 
same example to demonstrate the action of D7 on a partition of [6]. Consider 
the partition of {1, 2, 3, 4, 5, 6} with the following grammatical labeling 

{1,3,6}{2,5}{4} 

&6 ^5 b± a . 

If we apply the substitution rule a — >• abj of G7 to the letter a, then we get 
a partition with a consistent labeling 

{1,3,6}{2,5}{4}{7} 

bg 65 64 67 a . 

If we apply the substitution rule 65 — > bj of G7 to the letter 65, then we get 
the following partition with a consistent grammatical labeling 

{1,3,6}{2,5,7}{4} 

be &7 64 a . 

Theorem 5.3 For n>l, the multivariate polynomial S n (a,b) is stable. 

Proof. From the grammatical labelings, we see that S n (a,b) is multiaffine. 
Note that S n (a, b) is multiaffine in a, b\, 62, ■ • • , b n with every term containing 
a as a factor. Since 

G n+ i = {a -> ab n+ i,bi 6 n+ i, 1 < i < n}, 

for each multiaffine monomial of S n (a, b) which is of the form ah, we have 

D n+1 (ah) = ab n+1 h + aD n+ i(h). 

It follows that 

S n+ i(a,b) = D n+ i(S n (a,b)) 

n 

= b n+1 S n {a, b) + b n+1 ^2 (S„(a, b)). 

i=l 

Define 

T n+1 =b n+1 (l + Y^d/db^j . 

Then we have S n +i(a, b) = T n+ i(S n (a, b)). 

We proceed to prove the stability of S n (a, b) by induction on n. Since a is 
stable, we only need to show that the linear operator T n+ \ preserves stability 
of multiaffine polynomials. 

Let 

71 

F= (a + w) Y[(bi + Vi). 
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Then wc have 

71 

T n+1 (F) = bn+iF + bn+^d/db^F) 



1=1 



n 

= b n+1 F + 6„ +1 F^ — 



. bi + v t 

z=l 



bn+lF f 1 ' 



r-f fo » + w i / 

To prove that T n+ i(F) is stable, we assume that a, w, 6, and Vi have positive 
imaginary parts for all 1 < i < n + 1. Consequently, 

e = i ■ 



6j + «t 

is nonzero since it has negative imaginary part. Moreover, each factor of F has 
positive imaginary part, and so does 6 n +i. This yields that F and 6 n +i do not 
vanish. It follows that 

T„+i(F) =b n+1 F^0. 

Hence T n+ i(F) is stable. In view of Lemma we see that T n +x(S n (a,b)) is 
stable. This completes the proof. I 

It is worth mentioning that we use the operator T n+ \ instead of D n+ \ in 
the above proof because the operator D n+ \ docs not satisfy the condition in 
Lemma [Q1 Take D 2 as an example. It can be seen that D2((a + w)(bi + u)) is 
not stable. Note that 

D 2 ((a + w)(h + u)) = 62(0(61 + u + 1) + w). 

Let a = ^,61 = I — 1, u = I — 1 and w = i. But we have D 2 ((a + vS)(b\ +u)) 
This implies that i?2((a + w)(6i + u)) is not stable. 
Next we prove the stability of -B„(x, y, z, u, v). 

Theorem 5.4 Forn > 1, £/ie multivariate polynomial B„(x, y, z, u, v) is stable. 

Proof. Let /„ = D n D n -\ . . . Di(xq). From the grammatical labelings, it can be 
seen that /„ is multiaffine. We proceed to prove the stability of Di n Dm-\ ■ ■ ■ Di (xo) 
by induction on n. The stability of xo is evident. 

Wc now assume that fi n -i is stable. Let us consider the actions of Dm-x 
and Z?2n ■ By using the argument in the proof of the stability of C n (x, y , z) , it can 
be shown that the operator _D 2 „_i preserves stability of multiaffine polynomials. 
This leads to the stability of fi n -\ since jm-\ = D2n-i{f2n-2) ■ 

Recall that 

G 2n = {xi,yi,Zi,Ui,Vi ->■ x n y n z n ,0 ^ i < n; 

^ x n z n u n7 v n y x n y n z n }. 
Let B denote the following alphabet 

{xi,Vi,Zi, u h Vi,0 < i < n} U {«„}. 
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Since /2n-i is multiaffine and each term in f2 n -i contains u n , we may write a 
monomial of fm-\ in the form u n h. Then we have 

D2n(u n h) = (x n z n u n )h + D 2n (h). 

Thus, 

fin = ^2n(/2n-l) 

— X n Z n j2n — 1 t X n y n Z n ^2 d/d w (f 2n -i). 

weB 

Hence we may write f 2n as T(/2n-i), where T is a linear operator as given by 

w£B 

It remains to show that T preserves the stability of multiaffine polynomials. 

Let 

F = (u„ + r Un ) (w + r w ). 

wEB 



By Lemma I5T] it suffices to verify the stability of the following polynomial 
T (F) = x n z n F + x n y n z n F V" — — 

^ — ^ in -4- 1 



wEB 



/l ^ 1 \ 

v n y„z n F h y — ; 



Suppose that all the variables Xi,yi, Zi,Ui,Vi,r Xi ,r yi ,r Zi ,r Ui and r Vi have posi- 
tive imaginary parts for < i < n. Then 

has negative imaginary part, and so it is nonzero. Meanwhile, every factor of 
F is nonzero since its imaginary part is positive. Note that under the above 
assumption, x n , y n and z n have positive imaginary parts, and hence they are 
nonzero. Consequently, T (F) = x n y n z n F^ does not vanish. This leads to the 
stability of T(F). 

In light of Lemma 15.11 we deduce that / 2n is stable. This completes the 
proof. I 

The proof of the stability of C n (x, y, z) applies to the stability of T„(x, y, z). 
The details are omitted. 



Theorem 5.5 For n 1, the multivariate polynomial T„(x, y,z) is stable. 



Multivariate stable polynomials can be reduced to real-rooted univariate 
polynomials by diagonalization and specialization, see Wagner [TS]. More pre- 
cisely, if / £ R[zi,Z2, -.-tZn] is stable, then/(zi, . . . , z n )\ Zi=Z] and/Oi, . . . 

: Zn j | zi—a 

are also stable, where 1 ^ i ^ j ^ n and o£t. For example, setting a = 1 and 
&i = &2 = ■ • • = b n = x in S n (a,b) leads to the real-rootedness of S n (x), see 
Harper [12]. 
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For the multivariate stable polynomials B n (x, y, z, u, v), applying the diag- 
onalization yi — Vi ~ x and the specialization xi = Ui — Zi = 1 for all < i < n. 
we are led to Theorem II .31 

Let M (n, k) denote the number of Legendre-Stirling permutations of M n 
with k barred descents. By setting Xi = u, = j/j = Zi = 1, and Vi — x for all 
< i < n, we obtain the real-rootedness of the generating function of M(n, k). 

Corollary 5.6 For n ^ 1, t/ie polynomial 

n 

M n (x) =^2M(n,k)x k 
fe=i 

/ias only real roots. 

For the multivariate stable polynomials T„(x, y,z), by setting x\ = Zi = 1 
and y,; = y for all ^ i ^5 n, we are led to the real-rootedness of T n (y), which 
is equivalent to the real-rootedness of C n {x). 
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